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The stationary nonlinear Schrodinger equation (or Gross-Pitaevskii equation) for one-dimensional 
potential scattering is studied. The nonlinear transmission function shows a distorted profile, which 
differs from the Lorentzian one found in the linear case. This nonlinear profile function is analyzed 
and related to Siegert type complex resonances. It is shown, that the characteristic nonlinear profile 
function can be conveniently described in terms of skeleton functions depending on a few instructive 
parameters. These skeleton functions also determine the decay behavior of the underlying resonance 
state. Furthermore we extend the Siegert method for calculating resonances, which provides a 
convenient recipe for calculating nonlinear resonances. Applications to a double Gaussian barrier 
and a square well potential illustrate our analysis. 
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I. INTRODUCTION 

Transport properties of Bosc-Einstcin condensates 
(BECs) are of considerable current interest, both experi- 
mentally and theoretically. Especially atom-chip exper- 
iments are well-suited to study the influence of inter- 
atomic interaction on transport properties of BECs in 
wavegiudes since different waveguide geometries can eas- 
ily be realized [l], [3, H, EI- An alternative method was 
implemented in a recent experiment [f| where a BEC was 
created in an optomagnetic trap and outcouplcd into an 
optical waveguide. 

A convenient theoretical approach is based on 
the Gross-Pitaevskii equation (GPE) or nonlinear 
Schrodinger equation (NLSE) 
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which describes the dynamics in a mean-field approxima- 
tion at low temperatures 0, 0, H, @ ■ Another important 
application of the NLSE is the propagation of electro- 
magnetic waves in nonlinear media (see, e.g., [Io| . ch. 
8). The ansatz ijj(x,t) = exp(— i^Lt/h) i/j(x) reduces ([1]) 
to the corresponding time-independent NLSE 



rare cases were the one-dimensional NLSE 
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can be solved analytically [3 US, HUES- As an example, 
Fig. Q] shows the nonlinear transmission coefficient \T\ 2 
as a function of the chemical potential for the square well 
potential considered in [T3 |. which is discussed in more 
detail in section ITVl 
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FIG. 1: Transmission coefficient obtained from solving the 
stationary NLSE Q for a square well potential with a = 20, 
) +.g|^(x)| 2 ) ^(x) =A*V(x) (2) Vb = -50 for g = +1 (see Sec.|IVB|and 0). 



with the chemical potential fi. 

Various interesting phenomena have been reported 
originating from the nonlincarity of Eq. |T]), as for in- 
stance a bistability of the barrier transmission probabil- 
ity [ll|, E2, EH, EI| • A paradigmatic model in this context 
is the transmission through a one-dimensional rectangu- 
lar barrier or across a square well potential, one of the 
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One observes a clear structured behavior: the well- 
known Lorentz profiles determined by the complex- 
valued resonances, well understood for linear transmis- 
sion, are distorted. The curves bend to the right (to 
the left for attractive nonlincarity g < 0) and arc multi- 
valued in certain regions. A time-dependent numerical 
analysis shows that the lowest branch of the transmission 
coefficient is populated if an initially empty waveguide is 
slowly filled with condensate with a fixed chemical po- 
tential fi [ll], EH whereas the highest branch can be pop- 
ulated if the chemical potential is adiabatically increased 
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during the propagation process which is equivalent to 
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plying an additional weak time-dependent potential 
These and related aspects are discussed in detail in [18 1. 

It is the purpose of the present paper to determine the 
functional form of the transmission profile surrounding 
resonance peaks characterized by \T\ 2 = 1 for a general 
(symmetric) potential by relating them to Siegert type 
complex resonances. In particular we show that the line- 
shape of the resonance peaks is determined by the decay 
behavior of the underlying mctastable resonance state. 

The paper is organized as follows. In Sec. [TTJ, we de- 
velop a formula for the nonlinear Lorentz profiles which 
describes the transmission coefficient in the vicinity of a 
resonance in terms of skeleton functions which also de- 
termine the decay behavior of the resonance. In Sec. IIII1 
we present a convenient recipe for calculating nonlinear 
resonances and skeleton curves which we call the Siegert 
method. In Sec. IIV1 we use this method to demonstrate 
the validity of the nonlinear Lorentz profile from Sec. |TT] 
for two example potentials. In Appendix [X] we discuss 
the continuation of solutions of the NLSE to complex 
chemical potentials and in Appendix [B] we derive a use- 
ful formula for the decay coefficient of a symmetric finite 
range potential which we call the Siegert relation. 



II. NONLINEAR LORENTZ PROFILE 



A. The transmission problem 



coefficient is nontrivial. Here we review and slightly ex- 
tend an approach based on the time-dependent NLSE 
(see [O, [lj, [Hi). Following Paul et al. [3 we con- 
sider an experimental setup where matter waves from 
a large reservoir of condensed atoms at chemical poten- 
tial /i are injected into a one-dimensional waveguide in 
which the condensate can propagate (see Figure (2). In a 
time-dependent approach the system is described by the 
NLSE 
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where the source term f(t)e~ l,lt l h 5{x — xq) located at 
x — xo emits monochromatic matter waves at chemical 
potential fi and thus simulates the coupling to a reservoir. 
The barrier potential V(x) is assumed to be zero for x < 
x . 

In the following we assume a constant source strength 
/(£) = /o and look for stationary solutions ip{x,i) = 
ip{x) exp(—ijj,t/H) of Eq. ^ arriving at 



Hr/,(x) = - — iP"(x) + V(x)tP(x) 
2m 
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Application of the integral operator 
lim 
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leads to 
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where we have introduced the notation 

ip(x) ■■ 



tp-(x) x < xo 
ip+(x) x > x 
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Since we have V(x) = 0, x < Xo and there is no incoming 
current from x — —oo the solution in the region x < xo 
is given by the plane wave ip-(x) = ip°_ cxjp(-ik-x) with 

k- = ^/2m(/i — g\ip°_\ 2 )/h and ipL(x) = —ik-ip-(x). In- 
serting this into Eq. ((7]) leads to 



FIG. 2: (Color online) At xo condensate with chemical po- 
tential /i is injected into the waveguide from a reservoir. The 
reservoir emits a plane matter wave in both directions into the 
guide so that an incoming beam with current ji is partially 
reflected (current j re f and partially transmitted (current jt) 
at the barrier potential V(x). 

In the case of the NLSE the superposition principle 
is not valid. Therefore the definition of a transmission 
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or, taking into account the continuity of the wavefunction 



at x 

r, 2 



Xo, to 
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with fc + = y/2m(n - g\^ + (x )\ 2 )/h. Eq. (fTO]) relates 
the wavefunction ip+ in the region x > xo to the source 
strength fo- In order to relate the source strength with 
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the incoming condensate current we consider the special 
case where V(x) = everywhere, i.e. without a barrier. 
Then the wave function in the region x > xq is given 
by a plane wave i/j+{x) = Aexp(ikA(x — xO)) with fc^ = 
\/2m(/i — g\A\ 2 )/h and ip' + (x) = —ikA^p+{x). From the 
continuity of the wave function ip_(xo) = ip + (xo) =iwc 



get k- = k + = kA and together with Eq. (|T0|) we obtain 
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For a given source strength /o, Eq. (fTTj) can have up to 
two different solutions for A. In the following we only 
consider the solution corresponding to the limit of weak 
interaction. The incoming current emitted by the source 
is given by j in = ^k A \A\ 2 = \fo\\A\/h. InsertingEq. JllD 
into Eq. (jTUJ) we obtain the relation 



2ik A A = (ip' + + ik + if; + ) = 



(12) 



connecting the condensate wavefunction with the incom- 
ing current. We define the transmission coefficient as 



\T\ 
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where the current jt transmitted through the barrier is 
obtained by evaluating the current operator 



Jt 
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anywhere in the region x > xq, and j\ n is the current 
in absence of the barrier. In the noninteracting limit 
g — ► this definition coincides with the usual defini- 
tion of the transmission coefficient known from the linear 
Schrodinger theory. It has the advantage of being appli- 
cable in the time-dependent case. An alternative way 
to extend the concept of transmission to the interacting 
case is discussed in jl8l ]. 



B. Resonance lineshape 

In the following we will derive a formula for the trans- 
mission coefficient \T\ 2 of a symmetric potential well 
V(x) — V{— x) with the finite range a (i.e. V(x) = 
if |x| > a) in the vicinity of a resonance in dependence 
of the chemical potential [i of the incoming condensate 
current. Our approach is based upon a generalization of 
Siegert's derivation of the dispersion formula for nuclear 
reactions and we closely follow the arguments in [l9| ■ An 
alternative ansatz makes use of the Fcshbach formalism 

MM- 

We consider the situation where the condensate source 
is located at xq = —a. The solution in the down- 
stream region x > a is then given by a plane wave 
tp(x) = C exp(ikcx) with kc = \/2m(/j — g\C\ 2 )/h so 
that the transmitted current is given by jt = —kc\C\ 2 . 



At x = —a the wave function must satisfy Eq. (|12[) with 

Xq = —a. 

Thus the scattering wavefunction ij> in the interval 
[—a, a] is a solution of Eq. ([3]) with boundary conditions 

ip(a) = Ce ikca , i;'{a) = ik c ^{a) (15) 
2ik A A = ip'(-a) + ikip(-a) = (16) 



where k = y2m(ju — g\tp(— a)\ 2 )/h. The transmission 
coefficient given by \T\ 2 = j t /j ln = (k c /k A )\C/A\ 2 de- 
pends on the chemical potential /x and, due to the non- 
linear term in ([3]), also on the magnitude of the wave- 
function. If the incoming amplitude A is kept fixed, this 
dependence can be conveniently described by the magni- 
tude \C\ 2 of the outgoing amplitude. 

Now wc consider C/A as a function of fi. From (fT5"l) - 
([TBI we obtain 



C 2ikAip(a) exp(— i(k + kc)a) 
A ikip(-a) + ip'(-a) 
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Singularities of (|17|) occur for certain complex chemical 
potentials W s k where the denominator vanishes. These 
values of the chemical potential are defined by the eigen- 
value problem 
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in — a < x < a with the boundary conditions 



^ S k(a) = Ce^ kQ , VskW = ^sk^k(a) (19) 

iM-a) = Ce^ a + iS \ i>' sk {-a) = -ik s ^ sk (-a) (20) 



where k^ = y / 2m(W sk — g\C\ 2 )/h and 8 is some real 
valued phase. Because of the nonlinear term in (fT5)l . the 
complex energy VF s k = /x 8 k — jr s i;/2 with real ^ s k and r s k 
depends explicitly on \C\ 2 . The problem concerning the 
continuation of the solution to the domain of complex 
chemical potentials is discussed in Appendix [SJ Moti- 
vated by the analogy to a driven nonlinear oscillator we 
call the functions /iskdC*! 2 ) and r s k(|C| 2 ) skeleton curves 
and V'sk the skeleton wavefunction. 

From Eq. (fT5|) and its complex conjugate as well as 
the boundary conditions (HHJ), (|2H)) we derive the useful 
formula 
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In order to obtain C/A in the vicinity of the singularity 
Wsk we multiply ([3J by ^ s k and (fl"8|) by tp and subtract 
these equations. By integrating the resulting equation 



g^m 2 - IVskl 2 ) = 
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from x = —a to x — +& and using the boundary condi- 
tions we arrive at 
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iip sk (a)ip(a)(k sk - kc) 
+tp ak (-a)(ik s -kip(-a) + ip'(-a)) 

+g f VskV'dV'l 2 - |V>sk| 2 )dx 
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+(W S k -/•*)/ ip&ipdx = . 
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Thus we can write the denominator of C/A in (fl"7|) as 

ip&ipdx 
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we obtain 



2m W sk — /i 
ft 2 fc sk + k c 

2m W sk - M -g(|g| 2 -|^(-q)| 2 ) 

ft 2 fc sk + k 



(24) 

(25) 
(26) 



ikip(-a) + ip'{-a) = 
W Bk - (j, 
(ft 2 /2m)V- sk (-a) 
2mg 



tp sk ^dx 



(27) 



- \^r)dx 



h 2 tp sk (-a) J_ a 

W sk - [i ip sk (a)il>(a) 

% 

(ft 2 /2m)f/' sk (-a) fc sk + k c 



.W f 



sk 



^-g(\C\ 2 -\i^a)\ 2 )^ sk {-a)^(-a) 



(ft 2 /2m)^ sk (-a) k sk + k 

Assuming that the eigenvalue W sk is not degenerate, we 
have in the limit fx — > W sk : %j) — > ip sk , \ip(— a)| 2 — > |C| 2 
and fc, fcc — > fc sk so that 



ikip(-a) + ip'(—a) 



W sk - /Lt 



ip sk dx + i- 
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For the numerator of C/A we get 2ikAip(a) exp(— i(fc -I 
kc)a) — > 2ifc^ sk (a) exp(— 2ifc sk a). Thus C/A becomes 
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For sufficiently small values of r sk . we can multiply ip sk 
by a suitable constant £ of magnitude 1 which makes CV'sk 
real (up to terms of order T sk ) in the regions of slowly 
varying phase which give the main contribution to the 
integral /" \tp sk \ 2 dx. In this limit we can thus use the 
approximation 



CV2k<fc« I \^ sk \ z dx. (30) 



We furthermore define the phase factor i5 sk by CV'sk(a) = 
| , sk (a)| exp(ifc sk a + iS sk /2). Using these definitions and 
V'sk(— a) = ip sk (+a) exp(iS) Eq. ([29)) can be written as 
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or, using ([2T 
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■-p 2m cxp(22/c s ka+i5 s k) (l+exp(i5)) 



As r sk tends towards zero, fc sk becomes real and the 
last term in the denominator of (|32p becomes negligible, 
so that we have in this limit 
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fcsk a* - Msk + «r sk /2 
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Thus the transmission coefficient in the vicinity of a res- 
onance is given by \T\ 2 = (kc Aa)|C '/ A\ 2 
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r s 2 k /4 



fc sk in - m*) 2 + r s 2 k /4 



(34) 



with fc sk f=a Zee 1 . This result formally resembles the 
Lorentz or Breit-Wigner form that occurs in the respec- 
tive linear theory. However the chemical potential /K sk 
and the width T sk depend implicitly on \T\ 2 . This depen- 
dence disappears in the linear limit g — > and we recover 
the usual Lorentz profile. Eq. (|3~4"|) can be inverted to the 
form 



, Tsk / kA 1 



(35) 



i. e. the skeleton chemical potential is the average of the 
two branches fi± and the skeleton r sk the (|T|-weightcd) 
width 



r sk = ( M+ (|T|)- M _(|T|)) 



in 



(k A /k sk ) - |T| 2 



(36) 
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Figure [¥] shows a typical nonlinear Lorentz curve of the 
type (|3~4"|) for the case of a repulsive nonlinearity. The 
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skeleton curve /i s t, indicated by the dashed line, appears 
as a kind of backbone structure of the nonlinear Lorcntz 
profile, justifying its name which is taken from the the- 
ory of classical driven nonlinear oscillators (see e.g. (2l| ) 
where resonance curves similar to (|34j) occur. 

The skeleton curves ^ s k and L s k can either be 
parametrized in terms of the amplitude |C| 2 or in terms 
of the number of particles N = J . |V'sk(^)| 2 ^ inside 
the potential well. It was shown (see e.g. [22|, [23j]) 
that in an adiabatic approximation the decay behavior 
of a resonance state is determined by the imaginary part 
r s k (N(t)) of the instantaneous chemical potential via 



d t N(t) 



r sk (N(t)) 



N(t). 



(37) 



Thus there is a close connection between the transmission 
lineshape and the decay behavior of the corresponding 
resonance state as it is known for the linear limit .9 = 
where the decay coefficient is constant and the lineshape 
is Lorcntzian. 



III. CALCULATING SKELETON CURVES 

As shown in Sec. [IT] the skeleton curves /Lt s k(|C| 2 ) and 
r sk (|C| 2 ) are obtained by solving the NLSE JTSJ) with 
the Siegert boundary conditions (|19p. ((20|) . It has been 
shown that the use of Siegert boundary conditions is 
equivalent to a complex rotation of the coordinates (see 
e.g. Hi!). Different procedures based on this principle, 
e.g. direct complex scaling or complex absorbing poten- 
tials, have been successfully applied to resonance states 
of the NLSE [22|, [H, Hfil • Here we present an alternative 
method which is numerically cheap, easy to implement 
and, though not quite as accurate as the complex scaling 
procedures, provides a convenient basis for approxima- 
tions. 

This method, which we call the Siegert method, is 
based upon neglecting the imaginary part — L s k/2 of the 
chemical potential and thus having only real values of 
fc s k = ^/27n(^ s k — g\C\ 2 )/h which is justified for not too 
large values of Y s ^. Since the boundary conditions (Tl^) 
and (f2"0"P can no longer be satisfied simultaneously for real 
values of fc s k, we replace the boundary conditions (|20[) by 
the less restricting condition 



ax 







x=0 



(38) 



which preserves the symmetry \ipsk(— x)\ 2 = |'0sk(^)| °f 
the skeleton wavefunctions. The new boundary value 
problem given by (TT8")) . fl"9")) and is solved using a 
shooting procedure where the NLSE JT5|) is integrated 
from x = +a to x = using a Runge-Kutta solver with 
starting conditions (|19p at x = a for a fixed value of C. 
By means of a bisection method, // s k is adapted until the 
condition (|3"8")l is satisfied. The imaginary part of the 



chemical potential can then be estimated by the Siegert 
relation 



sk 



h 2 k sk \cf 

Jo \Ak(x)\ 2 dx 



(39) 



where ±6 are the positions of the maxima of the sym- 
metric trapping potential (see Appendix |B|) . 

Before we use our simple method to compute skeleton 
curves in Sec. IIV1 we demonstrate its validity for the 
lowest resonance state of the standard test potential 



V(x) 



■ exp(— ax 2 ) 



(40) 



with a = 0.1 and b ~ 3.16 using units where h = 1 and 
m = 1 as we do for all numerical calulations in this paper. 
We choose a = 30 to be sufficiently large to ensure that 
the resonance wavefunction 4>sk(x) is well approximated 
by a plane wave in the area x > a . The amplitude 
C is chosen such that the wavefunction is normalized in 
the region |x| < b, i.e. J , \ip s k(x)\ 2 dx = 1. Note that 
because of the nonlinearity, ip s k is not proportional to C. 

In Table [J we compare our results for the lowest reso- 
nance of the potential (|40"]) with the results of direct com- 
plex scaling and the complex absorbing potential method 
|22l |. The agreement between the different methods is 
very good, especially if the interaction constant g is small 
which is also the case within our calculations of skeleton 
curves in the following section. Apart from being numer- 
ically cheap and easy to implement the Siegert method 
proposed here can provide analytical expressions for /z s k 
and L s k if the potential in consideration is simple enough 
(see subsection II V Bp . 



/'s 



0.4601 
0.7954 
1.0772 
1.3192 
1.5317 
1.7247 
1.9070 



0.4601 
0.7954 
1.0765 
1.3190 
1.5315 
1.7236 
1.9043 



0.4602 
0.7954 
1.0772 
1.3192 
1.5312 
1.7231 
1.9035 



r s /2 r cs /2 r C Ap/2 



9.63e-7 
1.81e-5 
1.56e-4 
8.11e-4 
2.79e-3 
6.78e-3 
1.27e-2 



9.35e-7 
1.82e-5 
1.55e-4 
8.05e-4 
2.76e-3 
6.65e-3 
1.24e-2 



9.62e-7 
1.80e-5 
1.56e-4 
8.05e-4 
2.75e-3 
6.63e-3 
1.23e-2 



TABLE I: Chemical potential and decay rates for the low- 
est quasibound state, calculated with the Siegert method 
(S), complex scaling (CS) and complex absorbing potentials 
(CAP). 



IV. APPLICATIONS 

In this section we compare exact results for transmis- 
sion peaks with the predictions of the nonlinear Lorentz 
profile for two different model potentials. To this end 
we calculate the skeleton curves /x s k and r s k, which can 
either be parametrized in terms of the amplitude |C| 2 or 
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in terms of the number of particles N = J. |^> s k(a:)| 2 cfcc 
inside the potential well, by means of the Siegert method 
presented in the previous section. Furthermore we show 
that these skeleton curves are conveniently approximated 
by polynomials depending on a few instructive parame- 
ters. 



A. Example 1: The double— Gaussian barrier 



To demonstrate the validity of our model (|34[) we apply 
it to the potential 

V(x) = V [exp(-(x + bf/a 2 ) + cxp(-(x - b) 2 /a 2 )] 

(41) 

with the parameters Vb = 1, b = 14.7/2, a = 6/5 and 
a nonlincarity of g = 0.005. In [TTI ] the transmission 
coefficient of this potential in dependence of /i is calcu- 
lated for the case of an initially empty waveguide. The 
incoming amplitude \A\ 2 is connected with the incident 
current j; n (i.e. the current in absence of the barrier) via 
ji n — \A\ 2 y / 2(fi — g\A\ 2 )/m. In the following we will as- 
sume \A\ 2 = 1 in all numerical calculations. 

Using the method described in Sec. IHII we numerically 
calculate the skeleton curves /^ s k(|C| 2 ), r s k(|C| 2 ) with 
\C\ 2 < \A\ 2 . For \C\ 2 = \A\ 2 we have |T| 2 = 1 (see 
Sec. HI1)- We call the quantities hr := n s k(\C\ 2 = \A\ 2 ), 
r n ■■= T sk (\C\ 2 = \A\ 2 ) and ^ H (x) := ^ sk (x; \C\ 2 = \A\ 2 ) 
the resonance chemical potential, resonance width and 
resonance wavefunction respectively. In the limit |C| 2 — > 
the influence of the nonlinear term in the NLSE ([18]) can 
be neglected so that /i s k(|C| 2 — ► 0) = fx n and r s k(|C| 2 — > 
0) = r„ where \x n and T n are the respective quantities of 
the linear problem with g = 0. 




FIG. 3: (Color online) The numerically calculated curves 
Afek(AQ, r s k(iV) and |C| 2 (iV) (solid blue) and the approxima- 
tion described by (|42p - (|46[) (dashed red) for the potential 
(|41[) as well as the Taylor approximation (|47[) (dashed dot- 
ted black curve) are almost indistinguishable on the scale of 
drawing. 




0.66 0.68 0.7 



0.72 0.74 0.76 0.78 



FIG. 4: (Color online) Comparison between the nonlinear 
Lorentz curve (|34[) (solid red curve) and the transmission co- 
efficient obtained from solving the stationary NLSE (dotted 
blue curve) for the double Gaussian potential (|41[) . Dashed 
blue: Skeleton curve /i s k(|T| 2 ). 



Now we will show that over a wide range of parameters 
the skeleton curves /i s k(A r ) and r s k(iV) are well approx- 
imated by simple elementary functions and that the five 
quantities r„, fiR, Tr and |^4| 2 provide all the nec- 
essary information. 

Since the shift in the chemical potential is caused by 
the term | "0sk | 2 in the NLSE (Tig]) , we assume this shift to 
be approximately proportional to the number of particles 
inside the potential well, that is 



Nr 



N 



(42) 



where Nr = f , \tpR(x)\ 2 dx is the norm inside the well in 
the case of resonance. Next we represent the amplitude 
| C | 2 = |i/)(a)| 2 as a function of TV by a Taylor series which 
we truncate after the quadratic term, 



|Cf« \T\ 2 \A\ 2 ^ xiN + X2N 2 



(43) 



If there are no particles (N = 0), the transmitted am- 
plitude \C\ 2 is zero so that there is no constant term in 
P3| . Inserting flU]) and (g3J) into J33J) we obtain 



r sk (A0 

2y/2h 



2h 2 k sk \Cf 



N 



(44) 



V^(N)-g(xiN + X2N 2 )( X i + X2N). 



From Eq. (|44|) we obtain 



r sk (iv = 0) = r„ = 2V2hy/ fin/mxi 

Nr = 2HV2^r~ 



(45) 

g\A\2\A\ 2 /(V^r R ) (46) 



so 



that 

Xi = r « 
Inverting 



the coefficients 
m/(2h 



Xi 



and 



X2 are given by 
Xi I Nr. 

Eq. (HU leads to N± = -Xi/(2 X2 ) ± 
Thus we can compute fi^ and 



2fi n ) and X 2 = \A\ 2 /N R 



VxWx¥+WJx2. 

r s k as a function of |T| 2 
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It is often useful to approximate Eq. (|4"4"|) by a second 
order Taylor polynomial 



r sk = r„(i + m N + V2 N 2 ) 



(47) 



where 771 = X2/X1 + (M-R ~ Pn)/{2p„N R ) ~ gXi/Pn 
and r? 2 = X2(m« - p, n ) / {2^i n N rXi) - gXi/Pn ~ 

[(PR - Pn)/(p n N R ) - gXl/Pn] ■ 

Figure [3] reveals an excellent agreement between the 
numerically calculated skeleton curves and the approxi- 
mation described by (02) - pj]) . 

Figure Q] compares a peak of the transmission coeffi- 
cient with the resonance model ([34]) . While the qualita- 
tive features such as the bending of the curve are well 
reproduced, the resonance model (J34J) slightly overesti- 
mates the width of the resonance curve. 



B. Example 2: The square well 

For illustrative purposes we now apply the result 
(|34|) to a simple analytically solvable toy model system 
which has a similar transmission behavior as the double- 
Gaussian barrier considered in the previous section. In 
addition it shows resonance peaks originating from the 
bound states of the corresponding linear (g = 0) system 
which have been destabilized due to repulsive (<? > 0) 
interaction and have thus undergone a transition from 
bound to resonance state (3]. We consider the finite 
square well potential with vanishing interaction outside 
the potential well, where the wavefunction ip(x) must sat- 
isfy 



Am (far 



and 



2m dx 2 



+ g\^{x)\ z + V o -p)^(x)=0, \x\<a. 



(49) 

This model with vanishing interaction outside the po- 
tential well was introduced in [l4[ in order to discuss 
the scattering process in terms of ingoing and outgoing 
waves and thus enabling an analytical treatment. In prin- 
ciple the interaction outside the potential well can be 
eliminated by means of a magnetic Fcshbach resonance 
(see e. g. [27|) or by a larger transversal extension a± 
of the waveguide in this region since the effective one- 
dimensional interaction strength is proportional to 1 / a\ 
(see e. g. (f|). Alternatively, instead of neglecting the 
interaction outside the potential well one might add ad- 
ditional repulsive barriers at x = ±a without affecting 
much the qualitative behavior of the system with the 
disadvantage of making the analytical treatment more 
complicated. 

In [l4[ the transmission coefficient is calculated analyt- 
ically and it is shown that the respective states satisfying 



the NLSE ([H, gg]) with the boundary conditions (HD 
and (skeleton states) have the chemical potential 



Msk = V + -g\C\ 2 + 



h 2 K 2 {p)n 2 
2ma 2 



1 



- P 
2p 



(50) 



where n is an integer number and K{p) is the complete 
elliptic integral of the first kind (see e.g. [II])- The upper 
and lower alternative correspond to C| 2 ( Vb + <?| C| 2 ) > 
and < g\C\ 2 < \Vo\ respectively. The parameter < 
p < 1 is determined by 



g \A\ 2 [V + g\C\ 2 ]^=K\p) 



P 



(51) 



h 4 n 4 " ^' y p (p-i) 

and the norm of the wavefunction inside the well reads 



N = 2a\C\ 2 + 



2h 2 K(p)n 2 
gma 



K(p) - E(p) 



(52) 



Since g = in the region |x| > a we have fc s k = 
v / 2m/x s k/7i so that the decay width is given by 



2h^2jI R -\C\< 
mN 



(53) 



where /i s k and N are given in ([50)) and (T52^) . As in 
Sec. IIV Al the skeleton curves can be approximated by 
a Taylor polynomial. Applying our model to resonances 
of the square well potential (|48|) . (|49|) it turns out that 
it is sufficient to truncate the Taylor polynomials in ([4^]) 
and (gUl after the linear term. This leads to N/N R = \T\ 2 
and thus 



r sk 



Pn 

r„ 



{PR ~ Pn)\T\ 2 .. 

(r fl -r n )|r| 2 



(54) 
(55) 



Figures [5] and [6] show the transmission probability 
|T| 2 (/i) in the vicinity of a resonance, the exact solution 
and the resonance approximation introduced in Sec. IIIII 
for a deep square well with Vb = — 50 and a = 20. 
For both repulsive and attractive interaction, where the 
curves bend to the right or left, respectively, a good agree- 
ment between the nonlinear Lorentz curve (|34p in first 
order approximation (|54[) . (|55[) and the respective reso- 
nance peak (see 14[) is observed. In particular, Fig. [6] 
shows that the nonlinear Lorentz curve (|54"|) is also able 
to describe the unusually shaped peaks surrounding res- 
onances which correspond to bound states in the linear 
limit g — > (see [HI)- The deviations are due to the fact 
that only a single resonance is included in the present 
approximation. 



C. Decay behavior 

As discussed in section [TT] the decay behavior of the 
resonance state is described by 



d t N 



r sk (AQ 
h 



N . 



(56) 
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FIG. 5: (Color online) Nonlinear Lorentz curve (|34[) (solid red 
curve) and transmission coefficient obtained from solving the 
stationary NLSE (dotted black curve) for a square well poten- 
tial with a = 20, Vb = — 50 for the resonance with quantum 
number n = 129 for g — +1 (upper panel) and g — —1 (lower 
panel) . 



For the simple situation when the skeleton curves are ap- 
proximately linear in \T\ 2 , as in (|54|) . (|55|) . an analytical 
expression for the time dependence can be derived. With 
r sk (iV) = r„ + (T R - T n )(N/N R ) we obtain 



N 



-r„ 



N 



N R "Nr 

Separation of variables yields 

hdy 



(r R - r„) 



N 



dt 



T„y + (T R - T n )y 2 
with y = N/Nr. This can be integrated to give 

T n N R 



N(t) = 



r H (e r «(*-*o)A - 1) + r. 



(57) 



(58) 



(59) 



In the limit g —> this reduces to the linear decay 
behavior N(t) = Nr exp[— T n (t — to)/h}. In the limit of 
long times t S> io the system shows a linear decay N(t) — > 

[r„/r fl j7v K exp[-r„(t - t )/h] - exp[-r n (t - t Q )/h] as 

well. 

Figure [7] shows the decay according to formula (159")) 
and compares it with the numerical solution of Eq. (I56p 



0.8 



0.6 



o.4<<: 



0.2 



0.2 



0.4 



FIG. 6: (Color online) Nonlinear Lorentz curve (|34|l (solid red 
curve) and transmission coefficient obtained from solving the 
stationary NLSE (dotted black curve) for a square well poten- 
tial with a = 20, Vb = —50 for the resonance with quantum 
number n = 127, which had been a bound state in the linear 
case g = 0, for g = +1. 



and the linear decay in a scmilogarithmic plot. Formula 
(|5"9"|) agrees well with the numerical solution of Eq. ([Sri)) . 
In the limit of long times both curves are parallel to the 
linear decay curve so that the system adopts a linear 
decay behavior as predicted above. 
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20 



30 



t 



FIG. 7: (Color online) Decay according to formula (|59(l (solid 
red line), the numerical solution of Eq. ([56} (dashed blue line) 
and the linear decay behavior (dashed-dotted black line) for 
the parameters a = 20, Vb = — 50 and g — +2 . 



V. CONCLUSION 

We have presented an analysis of the nonlinear reso- 
nances found for transmission of a BEC through a one- 
dimensional potential barrier in a mean-field GPE de- 
scription. The Sicgcrt method for determination of res- 
onances is generalized to the nonlinear case providing a 
convenient recipe for the computation of nonlinear reso- 
nances. 

Based on this Siegert method, we developed a formula 
for the nonlinear Lorentz profiles which can be described 
in terms of skeleton functions depending on a few instruc- 
tive parameters. The skeleton curves also determine the 
decay behavior of the underlying resonance state thus 
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relating the transmission lineshape to the resonance life- 
time. 

Applications to a double Gaussian barrier and a square 
well potential illustrate and support our analysis. Finally, 
for a simple model an analytical expression for the decay 
behavior could be derived. We are therefore hopeful that 
the theoretical ideas presented may be useful in future 
work on nonlinear resonances. 
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APPENDIX A: ANALYTICAL CONTINUATION 



Since the NLSE 



2m 



tjj"(x) + V(x)ip(x) + g\^{x)\ 2 i){x) = fiij;(x) (Al) 



explicitly contains the squared magnitude |^>(a;)| 2 of the 
wavefunction it is not analytical and therefore the ana- 
lytical continuation of its solutions for complex values of 
/i is nontrivial. Following the arguments given in [29[ 
we decompose the solution of the stationary GPE as 
ip(x) = X(x)+iY(x) with real functions X(x) and Y(x). 
From the real and imaginary part of Eq. (|A1|) we get a 
system of two equations 



2m 



X"{x)+V(x)X{x)+g(X 2 (x)+Y 2 (x))X(x) 



/j,X(x) 
(A2) 



h 2 

2m y ' 



-V(x)Y(x)+g(X 2 (x)+Y 2 (x))Y(x) 



fiY(x) 
(A3) 

which are analytical. The solutions of this system of 
equations therefore have a straightforward continuation 
into the domain of complex chemical potentials. As 
an example, we consider the plane wave solution of the 
fr ee (V(x) = 0) GPE with ip(x) = Cexp(ik c x), k c = 
-\/2m(/i — g\C\ 2 )/h and C = \C\exp(i8). One can easily 
verify that its decomposition X(x) = \C\ cos(kcx + 6), 
Y(x) = \C\ sm(k c x + 6) satisfies the system (|A"2j) . (TA"3|) 
for all complex values of p. 



APPENDIX B: THE SIEGERT RELATION 

Derivation: In the following we will derive a formula for 
the decay coefficient of a resonance state of an arbitrary 
symmetric finite range potential. For the sake of gen- 
erality and for future applications we consider the cases 



of one, two and three dimensions simultaneously. For 
now we only consider resonances of the linear Schrodingcr 
equation (Eq. |T]) with g = 0). The applicability to the 
nonlinear case is discussed separately further below. 

Any solution ip(x,t) = \J p(x, t) exp[i0(x, i)] of the 
Schrodingcr equation (JTJ) with real functions 0(x, t) and 
p(x, t) > satisfies the continuity equation 



where 



d t p + divj = 



j = -pV0. 

m 



(Bl) 



(B2) 



Application of the Gauss theorem for vector fields leads 
to 



where 



dt.N = 



N 



p(x, t) d 1 



(B3) 



(B4) 



is the norm of the wavefunction within an D-dimcnsional 
volume V and A = d(V) is the directed surface of V. 
If the wavefunction is trapped inside the volume V the 
decay coefficient can be defined by the relation 



d t N = --N. 



Together with Eqs. (|B3]) and (|B4j) this leads to 



r = -h 



d t N 
N 



= h 



Jj-dA 



f v p(x,t) d^x 



(B5) 



(B6) 



Now we consider a radially symmetric potential 
V(x) = V(r) with finite range a, i.e. V(r) = if r > a, 
where r = |x|. We assume the potential V(r) to have 
a single maximum located at b < a. Assuming that the 
wavefunction varies slowly in time, we replace the time- 
dependent wavefunction ip(x, t) by the adiabatic reso- 
nance state - 0( x ) °f the stationary Schrodingcr equation 
( Eq. with g = 0). Due to symmetry, the wavefunc- 
tion in the area r < a can be written in polar coordinates 



^(x) = R(r)Y(fl) exp [i<f> r (r) + i<fa(n)] , r < a, (B7) 

with real functions R(r) and Y(£l) where f2 stands for 
the angle variables. The resonance wavefunctions of such 
a potential are obtained by applying purely outgoing 
(Siegcrt) boundary conditions, 



, a \ (D-l)/2 

1>(x) = R(a) (-) Y(Q) 



(B8) 



x exp [ik(r — a) + i4> r (a) + z</>n(f2)] , r > a, 



where k = Hc(y/2mfi/h). For narrow resonances where 
r/2 = — Im(^t) is small compared to e = He(fi) we can 
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make the approximation k ps y / 2mRc(/x) / ft. For D = 1 , 3 
(|B8|) is an exact solution, for D = 2 it only holds in 
the limit a — > oo (see below). This ansatz makes the 
wavefunction continuous at r = a. The continuity of the 
derivative implies the conditions 



la 



4>'M = k, 



(B9) 



where the prime denotes the partial derivative with re- 
spect to r. 

The resonance wavefunction shall be trapped in the 
region < r < b. Thus the volume V is a D-dimensional 
sphere with radius b so that d D x = r D ~ 1 drd£l and dA = 
e r b D ~ 1 dtt where e r is a unit vector in the radial direction. 
For the integral (|B3|) we need the scalar product e r • j = 
(h/m)R 2 (r)Y 2 (n)e r ■ V(0 r + fa). The volume integral 
(IB4I) becomes 



N = [ R 2 {r)r D ~ 1 dr f Y 2 (Q)dn. (BIO) 
Jo Jn 

For the surface integral (|B3|) we make the approximation 



d t N = 



[ 2-dAf* [ j-dA' (Bll) 

J A J A' 



where A' is the surface of the sphere with radius a which 
means that the reflection in the region b < r < a is 
neglected and (|B3|) becomes 



-9tNm [ ydA' = —a D - 1 R 2 (a)(f>'(a) / Y 2 (Q)dtt. 
J A' rn Jn 

(B12) 

By inserting Eqs. (|Bl2]) . (|B10|) and (jB9]) into (JB6]) we 
finally obtain the formula 



r = 



ft 2 fc i? 2 (a)a 



D-l 



J° R 2 {r)r D - 1 dr 



(B13) 



for the decay coefficient of a resonance with the chemical 
potential fi of the potential V(r) with the finite range a 
which is trapped inside the region < r < b. 

The formula (jBHf fo r D = 1 resembles Eq. (|2"T]1. 
In contrast to formula (|B13[) . the wavenumber fc s k in 
Eq. (|21[) is complex and the integration extends over the 
whole region r < a. Thus Eq. (|B13[) can be regarded 
as an approximation to Eq. (|21[) (respectively to its 
generalization to higher dimensions (see 0))- 



In the two— dimensional case the ansatz (|B8[) is only 
an approximation. As promised above we discuss this in 



more detail now by inserting the ansatz 

ip(r) = R{r) exp{i<j> r {r)) (B14) 
with real functions R(r) and <p(r) into the radial part 
d 2 ip{r) 1 dip(r) 2m/z 



+ -^fi>(r)=0 



(B15) 



dr 2 r dr h 2 
of the two-dimensional Schrodingcr equation. Separating 

real and imaginary parts we arrive at 

R»(r) , ^(r) 0;2(r) + 2^ = Oj 



R(r) rR(r) 



h 2 



2 J^ + ±y> r (r) + K(r) = 0. (B17) 



The choice R(r) = C/y/r and </> r (r) = fcr + 0o with real 
constants C, fc and 0o solves the lower equation. The 
remaining equation yields k 2 = Im^jfi 2 + r~ 2 /4. 
Thus Eq. (|B15p is approximately solved by 
ip(r) = C / y/r exp(ikr + 4>o) with k = ^/2m/x/S if 
r 2 > ft 2 /(8m/x) so that the length a in Eq. (|B8[) must 
be chosen accordingly. 

Applicability to the NLSE: In the case of the NLSE 
([2]), the ansatz JB7|), ([B8j) and thus formula (|B13|) are 
still valid in many cases where the wave equation can 
still be separated into a radial part and an angular part. 

For D = 1, inserting the ansatz (|B7|) . (|B8|) into the 
nonlinear term in © leads to ^^(x)! 2 = gR 2 {r)Y 2 (£l) = 
gR 2 (r) since Y(fl) = ±1. This means that the non- 
linear term only modifies the radial part of the wave 
equation whereas the angular part is not affected. For 
r > a, g|?/»(x)| 2 = gR 2 (a) is a constant term which only 
causes a shift in the chemical potential. If gR 2 (a) < 
Re(/i) this can be accounted for by replacing /.t — > 
/i — gR 2 (a) so that the wavevector is now given by 
k = Re( x /2m( f i- gR 2 {a))/h). 

If D = 2, we also have 1^(0) | = 1, so that the wave- 
function can still be separated into a radial part and an 
angular part in analogy to the one-dimensional case. For 
r > a, g\^p(x)\ 2 = gR 2 (a)a/r. Thus for \g\R 2 (a) < |pt| 
we can neglect the influence of the nonlinear term in the 
region r > a and the ansatz (|B8|) is still valid. 

For D = 3, we have |^(^)| 7^ 1 in general so that the 
wave equation can no longer be separated into a radial 
part and an angular part an thus the ansatz (|B7|) . (|B8|) 
fails. However, for the special case of s-wave solutions 
|F(S!)| = 1 still holds and in analogy to the case D = 2. 
formula (|B13|) is approximately valid if \g\R 2 (a) <C \n\. 
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